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Light can exert radiation pressure on any object it encounters and that resulting optical force 
can be used to manipulate particles. It is commonly assumed that light should move a particle 
forward and indeed an incident plane wave with a photon momentum hk can only push any 
particle, independent of its properties, in the direction of k. Here we demonstrate, using 
full-wave simulations, that an anomalous lateral force can be induced in a direction per- 
pendicular to that of the incident photon momentum if a chiral particle is placed above a 
substrate that does not break any left-right symmetry. Analytical theory shows that the 
lateral force emerges from the coupling between structural chirality (the handedness of the 
chiral particle) and the light reflected from the substrate surface. Such coupling induces a 
sideway force that pushes chiral particles with opposite handedness in opposite directions. 
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Electromagnetic (EM) waves carry linear momentum as each 
photon has a linear momentum of M in the direction of 
propagation. Circularly polarized light carries angular 
momentum due to the intrinsic spin angular momentum 
(SAM) of photons 1-6 . When light is scattered or absorbed by a 
particle, the transfer of momentum can cause the particle to 
move. Thus light can be used to manipulate particles 7-25 . Light 
will push a particle in the direction of light propagation (as 
illustrated in Fig. la) irrespective of the polarization of light and 
irrespective of the particle's own properties, even if it has chirality 
(Fig. lc,e), as long as we have a plane wave incidence (that is, a 
well-defined k). Let us now consider the configuration shown in 
Fig. lb, which shows a particle placed close to a surface made of 
an ordinary material (for example, Au or Si or silica). In this case, 
one might still expect the particle to be pushed in the direction of 
light propagation, as the surface does not break the left-right 
symmetry. We will show that this is true only if the particle is 
non-chiral. If the particle has chiral character, however, it will 
experience an additional lateral force in a direction that depends 
on its own chirality as shown in Fig. ld,f. This counter- intuitive 
force comes from a lateral radiation pressure and an optical spin 
density force that couple the chirality of the particle to both 
the lateral linear momentum and SAM generated by the 
scattered wave of the chiral particle. The time-averaged SAM 
densities, defined as <L e > =s 0 /(4(ot)(E x E*) and {Lm)^^/ 
(4cl)/)(H x H*), respectively, for the electric and magnetic parts, 
are associated with the polarization of light 26,27 . The lateral force 
could move particles with chirality of different signs in different 
directions as shown in Fig. ld,f. A good example of a chiral 
particle is the helix shown in Fig. lg. Interactions between such 
kind of chiral objects and EM waves have been extensively 
studied 28 and are shown to give rise to interesting phenomena 
such as polarization conversion 29-34 , photonic topological edge 
states 35 and negative refractive metamaterials 36-39 . We will show 
that the coupling of EM near field and structural chirality (the 
handedness of a chiral particle) will induce an anomalous lateral 
force that pushes the particle sideways rather than forward. 

We first show the full-wave numerical results of the optical 
forces acting on helical gold particles induced by a linearly 
polarized plane wave. We then consider a simpler system 
consisting of a model chiral sphere placed above a substrate. 
To reveal the underlying physics and trace the origin of this 
intriguing phenomenon, we further simplify the configuration by 
considering a dipolar chiral particle above a substrate, in which 
case the problem can be analytically addressed. 

Results 

Lateral optical force on a gold helix placed above a substrate. 

Consider a gold helix with inner radius r = 50 nm, outer radius 
.R^lSOnm and pitch P = 300nm, as shown in Fig. lg. We 
are interested in the optical force acting on such a particle 
induced by a linearly polarized plane wave of the form 
E inc = zEq exp(ikoX — icot). The particle is located above a sub- 
strate with a gap distance of 10 nm and its axis is along the x 
direction (see Fig. ld,f). The substrate has the dimensions of 
Ix wx t (see Fig. lb) and it can be metallic (for example, gold) or 
dielectric. 

For an isolated helical particle, the scattering force induced by 
the plane wave will push it in the direction of ko independent of 
the handedness of the helix (see Fig. lc,e). However, if the helix is 
placed above a substrate (as shown in Fig. ld,f), an additional 
lateral force (F y ) will emerge and push it sideways. The lateral 
force acting on a left-handed (LH) and a right-handed (RH) 
particle takes opposite signs. The optical force is calculated 
numerically by the Maxwell stress tensor method (see Methods). 



Figure 2a shows the lateral forces acting on an LH (blue lines) 
gold helix and an RH (red lines) one consisting of four pitches 
when they are placed above a gold (described by a Drude model, 
see Numerical simulation section) substrate (solid circles) and a 
dielectric substrate with e d = 2.5 (hollow triangles). The lateral 
forces are evident in a wide range of frequencies, where local 
resonances associated with the geometric shapes of the particles 
result in some oscillations. It is important to note that the lateral 
force takes opposite signs for the LH and RH helices. The 
existence of a lateral force in the case of a perfect electric 
conductor substrate is also examined and similar results are 
obtained (see Supplementary Fig. 1). 

Figure 2b shows the dependence of the lateral force on the 
number of pitches there are in the helix at /= 490 THz. The 
magnitudes of the lateral forces in the cases of a gold substrate 
(solid circles) increase with the pitch number, which indicates 
that the magnitude of the lateral force increases with the chirality. 




Figure 1 | Anomalous lateral force in the helix particle-substrate 
configuration, (a) A linearly polarized plane wave pushes a normal 
spherical particle forward. A standalone RH helix (c) and a standalone LH 
helix (e) will be pushed forward by radiation pressure. When the normal 
spherical particle is put near a substrate, (b) it still moves forward while the 
helical particles (d,f) are shifted in opposite directions by an anomalous 
lateral force. The arrows associated with the helical particles indicate the 
handedness. The substrate has the dimensions of / x w x t. 
(g) Dimensions of the helix particle. It has inner radius r = 50nm, outer 
radius R = 150nm and pitch P = 300nm. 
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Figure 2 | Numerically calculated lateral forces acting on the helical gold 
particles, (a) Lateral forces F y acting on the LH (blue lines) and RH (red 
lines) particles in the presence of metal (gold) and dielectric (e d = 2.5) 
substrates, (b) Lateral forces as a function of the number of pitches for the 
LH (blue lines) and RH (red lines) particles in the presence (circles) and 
absence of (triangles) the gold substrate, (c) Lateral force acting on 
the LH particle as a function of the gap distance between the particle and 
the gold substrate, (d) Lateral force as a function of the thickness of the 
dielectric substrate for the LH particle. In the numerical simulations, the 
frequency is set at f= 490 THz for case of the metallic substrate (b,c) 
and at f= 380 THz for the case of the dielectric substrate. 

In the absence of a substrate, there is a small residual lateral 
force (solid triangles) due to the end effect and this residual 
lateral force decreases as the pitch number increases (the residual 
lateral force is zero if the particle is symmetric on the jyz-plane). 
Figure 2c shows that the magnitude of the lateral force decreases 
when the gap distance between the particle and the gold substrate 
is increased, indicating that the force is due to the coupling 
between the particle and the substrate. Figure 2d shows the lateral 
force as a function of the thickness t of the dielectric substrate, 
where the frequency is set at /= 380 THz. The lateral force 
undergoes oscillations and we will see later that the oscillations 
are due to the Fabry-Perot resonances associated with the 
dielectric substrate. Note that we do not need to consider the 
effect of thickness in the case of the gold substrate due to the 
finite penetration depth of fields. 

The numerical results are rather counter- intuitive. The 
standalone helix (that is, without a substrate), whether it be LH 
or RH and assuming it is long enough so that the end effect can 
be ignored, scatters the same amount of light to the left ( -\-y) and 



to the right ( — y) and hence light cannot push it sideways. 
However, when the helix is placed above a substrate, it 
experiences a lateral force in the y direction. To see pictorially 
how this happens, we examine the magnetic field (H x ) patterns 
for both an RH (Fig. 3a) and an LH (Fig. 3b) helix in the case of 
the gold substrate at /= 420 THz. Figure 3a,b shows that the 
scattered field is asymmetrically distributed. The surface waves 
propagate predominately in the +y direction in the case of the 
RH particle while they propagate in the -y direction in the case of 
the LH particle. These field patterns show that the helix-above- 
substrate configuration produces asymmetric scattering and 
hence a lateral force, although individually the helix and the 
substrate produce symmetric scattering. 

Lateral optical force on a chiral sphere placed above a 
substrate. The numerical results on the gold helix demonstrates 
the existence of the lateral force. We now show that the effect can 
be observed even for a chiral sphere described by simple con- 
stitutive relations. Consider again the configuration shown in 
Fig. lb, where the spherical particle with radius r = 75 nm is now 
made of a model bi- isotropic chiral material described by the 
constitutive relations 40-44 : 



D 








E 


B 




- nc/c fi Y fi 0 _ 




H 



(1) 



where e r and ja r are the relative permittivity and permeability of 
the material, respectively; k is the chirality parameter and takes a 
real value here; e 0 > an d c are the permittivity, permeability and 
the speed of light in vacuum. The dispersion for a plane wave in 
such a medium is k± = (y/s Y ^ r ± k)co/c, corresponding to 
circularly polarized states of light. The gap distance between the 
sphere (e r = 2, fi r —l) and the surface of the gold substrate is 
10 nm. The particle is under the excitation of the same linearly 
polarized plane wave defined before. The wavelength is set at 
X — 600 nm. Figure 4a,b shows the time-averaged Poynting vector 
for the cases with and without the gold substrate, computed by 
the numerical solver COMSOL Multiphysics (http:// 
www.comsol.com/). For the isolated particle (Fig. 4a), the 
Poynting vector has a rotation pattern due to the chiral 
property of the sphere and the scattered energy is symmetric 
with respect to -y and +y. Such rotating Poynting vector pattern 
can also be found in the case of the gold helix (see Supplementary 
Fig. 2). However, in the presence of the gold substrate, the 
'handedness' of the Poynting vector in the near-field region 
results in an asymmetric pattern as shown in Fig. 4b and hence a 
lateral force. Figure 4d,e shows the flow of the time-averaged 
electric SAM[<L e > =e 0 /(4coi)(E x E*)]. Similar to the Poynting 
vector distribution, the computed SAM has a symmetric 
^-component in the isolated sphere case (Fig. 4d). However, in 
the presence of a gold substrate shown in Fig. 4e, the SAM tends 
to point in —y direction. We will show (through analytical 
theory) that the asymmetric distributions of Poynting vector and 
SAM are the consequences of particle-substrate coupling and are 
closely related to the origin of the lateral force. The blue line in 
Fig. 4g shows the numerically evaluated lateral force as a function 
of the chirality parameter k for the case with the gold substrate. 
The magnitude of the force increases with the magnitude of the 
chirality. If the chirality is positive, the force is negative, and vice 
versa. The force vanishes when the medium is non-chiral (k = 0). 
The chirality determines the direction of EM energy coupling and 
this property can be used to realize uni- directional excitation of 
surface waves 45-47 . 

The results shown in Fig. 4a-f also provide an intuitive 
understanding of the lateral force. From a symmetry point of 
view, the time-averaged Poynting vectors corresponding to the 
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Figure 3 | Magnetic field distribution on the gold substrate, (a) H x field pattern for the RH particle on the gold substrate, (b) H x field pattern for 
the LH particle on the gold substrate. Note the asymmetrical pattern and the different directions in which the RH and LH particles scatter EM fields. The 
plotted fields are on a plane just above the substrate and the frequency is set at f=420THz. 
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Figure 4 | Numerical results for a chiral sphere above a gold substrate. Time-averaged Poynting vectors (a) for an isolated chiral sphere, (b) for a 
chiral sphere above a gold substrate and (c) for a chiral sphere sandwiched symmetrically by two gold substrates. We set k = 1 for the sphere. Time- 
averaged electric spin density <L e > (d) for the isolated sphere, (e) for the sphere above a gold substrate and (f) for the sphere sandwiched by two gold 
substrates. The left-right asymmetry is obvious in panels (b) and (e). (g) Lateral force acting on the chiral sphere above a gold substrate (blue line) 
and sandwiched by two gold substrates (red line) as a function of the chirality parameter k. The blue line shows that the sign of the lateral force F y depends 
on k and F y = 0 if k = 0. The red line indicates the lateral force vanishes in the sandwiched case. The frequency is set at f=500THz. 



total fields for a standalone chiral/spiral particle under linearly 
polarized light excitation must satisfy <S(jy, z)} — — <S(— y, 
— z)y due to the rotational symmetry of the system. This can be 
seen by examining the numerical results in Fig. 4a. This implies 
that the total momentum flux scattered by the chiral object to the 
right (—y) and left (+y) should be the same for a standalone 
particle, which in turn implies that there cannot be a net force in 
the lateral (y) direction. However, if we break the z— > — z 



symmetry of the environment by putting a substrate underneath 
the particle, there is no symmetry requirement for | < S(y, z)}\ 
and | < S( — y, — z) ) | to be equal. This is indeed the case for the 
numerically calculated Poynting vector pattern as shown in 
Fig. 4b. A lateral force (F y ) can now exist as the total photon 
momentum scattered to the left and right are not balanced, and 
the magnitude of the lateral force depends on the details of near- 
field coupling. If we recover the symmetry of the environment by 
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adding another substrate above the sphere, the ^-components of 
the Poynting vectors (Fig. 4c) and the SAM (Fig. 4f) become 
symmetrical and the lateral force vanishes again as numerically 
verified in Fig. 4g (red line). 

The results here show that the lateral force emerges in a very 
general configuration of a chiral sphere above a substrate. The 
sign and magnitude of the lateral force are both directly related to 
the chirality of the particle. In the next section, we show with the 
help of an analytical theory that the lateral force is related to the 
chirality- generated lateral photon momentum. 

Force on a dipolar chiral particle placed above a substrate. The 

results shown up to now are full-wave simulations but in order to 
obtain an intuitive understanding of the origin of the lateral force, 
we consider the configuration of a dipolar chiral particle above a 
substrate in which case the lateral optical force can be analytically 
evaluated and we will show that it comes from the lateral 
radiation pressure and spin density force arising from reflection. 
The induced dipole moments of such a chiral particle can be 
expressed as follows 28 : 
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where p and m are the electric and magnetic moments, 
respectively; E and H are the fields acting on the particle. Here 
the polarizability of the particle is specified by parameters a ee , 
^mm and a em . We note that a em is related to the chirality 
parameter k of the material that the particle is made of and a em 
will changes sign if k changes sign and a em = 0 corresponds to a 
non-chiral particle. 

The optical force acting on a dipolar chiral particle in an 
EM field can be derived (see Supplementary Note 1) based on the 
expression of the force acting on a dipolar particle 48-50 , 
(F) = l/2Re[(VE*) -p + (VH*) •m-c^/(67r)(pxm*)], 
where E and H are the fields acting on the particle. The force 
expression can be written as follows: 

(S) 

(F) =Vt7 + (T^ -Im[a em ]Vx(S) + co- e Vx(Le> + c ( 7 m Vx(L m ) 
c 



+ wy e (L e )+coy m (L m )- 
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^Im[a ee a; m ]lm[ExH*], 



(3) 

where U= l/4(Re[a ee ]|E| 2 + Re[a mm ]|H| 2 — 2 Re[a em ]Im[H • E*]) 
is the term due to particle-field interaction; <S>=l/2 
Re[E x H*] is the time-averaged Poynting vector; <L e > and 
<(L m ) are the time-averaged spin densities defined before; a e — 
k 0 /s 0 Im[a ee ], G m = k 0 /[i 0 Im[a mm ] and a = a t + o m 

— c 2 ^/ (6n)(Re[(X ee (X^ m ] + OL em oil m ) are the cross-sections; y e = 

— 2cL>Im[a em ] + c/cQ/(37ieo)Re[a e ea* m ] and y m — — 2coIm[a em ] 
+ c/cQ/(37i/i 0 )Re[a mm aJ m ] also have the dimension of a cross- 
section. 

The first term in equation (3) corresponds to the gradient 
force. The second term represents the radiation pressure. The 
third term is a Vortex' force determined by the energy flow vortex 
around the particle and the optical activity (a em ). The fourth and 
fifth terms are scattering forces associated with the curl of the 
spin densities and are called curl-spin forces 51-55 . The sixth and 
seventh terms are referred to as spin density forces as they are 
directly related to the spin densities. 

For an isolated chiral particle acted upon by a linearly polarized 
plane wave, the gradient force and the Vortex' force vanish 
because the fields' amplitudes and the time-averaged Poynting 
vector are constants. The last term contributes nothing as 
Im[E x H*] = 0 for a plane wave. The spin density terms vanish 
as we are considering linear polarization. As a consequence, the 



expression of the force reduces to <F> = a < S > /c, which is just a 
forward-scattering force. 

In the presence of a substrate, the fields in equation (3) consist 
of both the incident plane wave and the reflected fields; that is, 
E = E inc + E ref , H = H inc + H ref , where the reflected fields can 
be expressed as follows: 
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Here G re f is the Green's function for reflection describing the 
effect of the substrate and it takes into account the reflected fields 
including the propagating and the evanescent components 42 . We 
now refer to equations (3-5) to reveal the origin of the lateral 
force in the presence of a substrate. In equation (3), no lateral 
component can come from the gradient force as the free energy 
does not change when we apply a virtual in-plane (xy-plane) 
displacement of the particle. For a lossless particle under long 
wavelength condition, its polarizabilities are dominated by their 
real parts and as a consequence, the strengths of the vortex force 
and the last term are much smaller than that of the spin density 
force. For example, the strengths of these two terms are, 
respectively, about 10 ~ 6 and 10 _ 3 smaller than that of the 
spin density force for a dipolar particle (e r = 2, /u r = 1, k= 1) of 
radius a — 30 nm placed at d = 60 nm above a semi- infinite gold 
substrate under the excitation of z-polarized plane wave 
(2 = 600nm). We can hence focus on the radiation pressure 
term and the terms related to spin densities. In the presence of a 
substrate, the lateral component of the Poynting vector becomes 
asymmetric (Fig. 4b) and hence contributes to a lateral radiation 
pressure. The magnetic spin density is small compared with the 
contribution of its electric counterpart because the non-magnetic 
substrate gives a relatively stronger response to the electric field. 
Hence the total spin density can be written as 
(L) = (L e ) + (L m ) n (L e ) « £o /(4«/)Im[E inc xE* ef ]. As the 
incident electric field E inc is along the z direction, this 
expression indicates that the direction of the spin density is 
mainly parallel to the surface of the substrate (as the numerical 
results shown in Fig. 4e confirm). Such a spin density can only 
contribute to a lateral force through the spin density terms but 
not through the curl-spin terms. 

As the spin density is mainly attributed to <L e ), the lateral 
component of the spin density force can be written as 
F y « coy e y • (L e ) ~ (y e s 0 /2)y • Im [E inc x E* ef ] . Since the incident 
electric field is polarized along the z direction, the lateral force 
requires an x- component reflected field. Such a component exists 
only when there is a substrate and it emerges from the 
interference between the particle and the substrate that enables 
the cross-coupling of different Cartesian components of the 
scattered fields. This can be understood by directly examining 
the expression of the reflected field. Consider the excitation of 
the chiral particle under the incident magnetic field H inc which 
is along y direction. According to equation (2), an electric 
dipole moment — /oc em H inc can be induced, which produces 
the reflected field E y e{ = ifi 0 tt e mW 2 G yy ei H[ nc in the presence of a 
substrate. Here G yy d is the corresponding tensor element of the 



Green's function for reflection. The reflected field E y e{ further 

^em-Eref 



induces a ^-direction magnetic dipole moment m y — — za e] 
and this magnetic dipole moment produces the required 
x-component electric field: E* e{ — ij2 0 a> 3 (xl m G yy f (V xG re f) xy H inc . 
Note also that (L e ) oc Im[Ei nc xE* ef ] oc a^and hence the spin 
direction remains unchanged when a em changes sign. Conse- 
quently, the sign of the lateral spin density force is determined by 
the sign of the coefficient y e as y e oc a em . 
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Consider a spherical dipolar chiral particle (a = 30nm) 
characterized by s r — 2.0 y /i r —l.O and a chirality parameter k 
and the particle is placed at d — 60 nm above a semi- infinite gold 
substrate. One can derive the polarizability tensor elements a ee , 
a mm and a em for the particle from its Mie scattering coefficients 
(see Supplementary Note 2). The lateral force contributed by the 
radiation pressure and spin density forces can then be analytically 
evaluated with the help of equation (3) (see Methods). Figure 5a 
shows the lateral force as a function of k at wavelength k — 600 
nm. It is clear that the lateral force versus chirality relationship is 
similar to those calculated using full-wave simulations for the 
bigger chiral sphere as shown in Fig. 4g. The analytical theory can 
also explain the oscillation phenomenon of the lateral force in the 
gold helix configuration (Fig. 2d). Let us now replace the gold 
substrate by a dielectric substrate (e d = 2.5 + 0.001/). We show in 
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Figure 5 | Analytical results for a dipolar chiral particle above a 
substrate, (a) Lateral force acting on a dipolar chiral particle ( a = 30nm, 
e r = 2.0) as a function of its chirality when the particle is located 60 nm 
above a semi-infinite gold substrate, (b) Lateral force acting on the chiral 
particle (k= -1) as a function of the thickness of a dielectric substrate 
(fid = 2.5 + 0.001/), showing oscillating behaviour, (c) Magnitudes of the 
reflection coefficients for an evanescent channel k// = 10k 0 (red line) and 
a propagating channel k// = 0.5k 0 (green dotted line and blue line). The 
wavelength is set to be A = 600nm. 



Fig. 5b the dependence of the lateral force on the thickness of the 
dielectric substrate according to the dipole theory. As expected, 
the force undergoes Fabry-Perot oscillations induced by the 
reflectance of the finite-thickness substrate. The dielectric 
substrate reflects the scattered field of the particle and causes 
constructive or destructive interference. A better understanding 
can be obtained by examining the reflection coefficient of the 
dielectric substrate described by the expression .R = 



[2/(C 2 -l)sinM]/[(C+l)e 



- ik z d 



(C- 



where C TE — 

kji/ijcoj, C™ = k z /(s d k 0z ) and k 2 z + kj / = e&n d 1%. Here k„ is 
the wave vector component parallel to the surface of the substrate. 
For a highly evanescent channel (fc//»fc 0 )» # TE = 0 and 
R™ « (e 2 d - 1)(1 - e 2 V)/[(e d - l) 2 - + 1)VV], where 
£™ is a monotonic function that increases with t and takes a 
limit value of (s d — l)/(s d -\- 1). The red line in Fig. 5c shows the 
magnitude \R™\ as a function of t for the case of /c//=10/c 0 , 
where we see that the reflection coefficient is a constant for large 
t. Propagating channels (k//<k 0 ), on the other hand, behave quite 
differently. Figure 5c shows the |i^ TE | (green dotted line) and 
|.R™| (blue solid line) for the case of ku — 0.5k Q . Clearly they 
undergo oscillations with a period that meets the condition 
k z tsX— 7i (here k = 600 nm), indicating the oscillations are due to 
the Fabry-Perot resonances. Combining the properties of both 
the evanescent and propagating channels, the final reflected field 
should oscillate and the resulting lateral force should too. 



Discussion 

While the lateral optical force is always there if the particle is 
chiral, it is clear that the lateral force is reasonably strong relative 
to the forward force only if the chirality is strong. The case with 
the helical gold particle, with numerical results shown in Fig. 2, is 
a good example in which the maximum value of the lateral force 
can reach about 0.4 times that of the forward scattering force if 
the gold helix is put above a gold substrate. We note that if the 
incident wavelength is much larger than the helix dimension, the 
lateral force would become very weak because the total scattering 
cross-section is small and in addition, the incident EM wave 
could not 'resolve' the geometric handedness detail of the helix 
which is crucial for inducing the chiral effect 56 . Although we have 
only shown that a linearly polarized plane wave can induce a 
lateral force on a chiral particle, the force also exists when the 
plane wave is circularly polarized (see Supplementary Fig. 3 for 
the case of gold helices on gold substrate). The key point here is 
the asymmetric coupling of a chiral particle with a substrate, 
which can also be induced if the incident light is the circularly 
polarized. In addition, the lateral force also can be induced when 
the incident plane wave propagates in a direction perpendicular 
to the axis of the gold helix (see Supplementary Fig. 4). 

In summary, we have numerically shown that the EM near- 
field coupling can induce a lateral optical force on a chiral particle 
near a substrate. The anomalous force can laterally push particles 
with opposite chirality in opposite directions. Using full-wave 
simulations, we established the existence of this counter- intuitive 
phenomenon for the helix-above-substrate configuration as well 
as the case of a chiral sphere with a simple constitutive 
relationship above a substrate. In the former case, the lateral 
force takes opposite signs for an LH and an RH helix. In the latter 
case, the lateral force changes sign when the chirality parameter 
changes sign from + k to — k. The asymmetric coupling between 
a chiral particle and a substrate breaks the left-right symmetry 
due to the special 'handedness' distribution of the Poynting vector 
(Fig. 4b). By analytically deriving and evaluating the lateral force 
acting on a dipolar chiral particle above a substrate, we show that 
such lateral force must exist even in the small particle regime, and 
the force is attributed to the lateral component of the radiation 
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pressure and 'spin density force' generated by the reflection field 
in presence of the substrate. The study reported here may find 
applications in the detection of chirality and in the separation of 
chiral molecules/enantiomers 57 using light induced forces. We 
note that the lateral force here must be accompanied by a torque 
induced on the chiral particles, the magnitude of which can be 
enhanced by a nearby surface. This may find applications in some 
opto -mechanical systems and as a light- driven 'rotor' or 
'motor' 58 . However, we note that thermodynamic fluctuation 
effects and hydrodynamic properties of the environment should 
also be considered in such applications. 

Methods 

The Maxwell stress tensor method. The optical forces acting on the gold helix 
and the chiral sphere are numerically evaluated using the Maxwell stress tensor 59 , 
which is denned as 



EE + c BB - 



1 



(E • E + rB • B)J 



(6) 



where I is the unit tensor. To calculate the forces, we first obtain the EM total field 
using a commercial finite-element-method package COMSOL Multiphysics (http:// 
www.comsol.com/). Then the tensor in equation (6) is integrated on a closed 
surface surrounding the particle. The forces shown in the main text are the time- 
averaged results. 

Semi-analytical evaluation of optical forces by dipole theory. When a spherical 
dipolar chiral particle is placed above a substrate, the fields of the induced dipoles 
will be reflected by the substrate and react on the particle. Hence, the induced 
dipoles can be expressed as 



Z(X em 
^mm 



Einc E r ef 
Hinc + H re f 



(7) 



For a small particle, the dipole polarizabilities in equation (7) can be related to the 
material parameters s r , fi r , k of the particle (see equation (1)) through Mie scat- 
tering theory (see Supplementary Note 2 and ref. 40). 

Combining equation (7) with equations (4) and (5), the reflected fields can be 
written as follows: 

E re f = £l - Einc + £2 " H inc , 
H re f — Ql * Ei nc + Ql ' Hi nc , 



where 



and 



Ci : 
C 2 - 

D 2 



(Ci -Dz + Ci) • (i-Cx-D^ 



(5i-e 2 +Di) • (f-Di 
(Ci -Di+C 2 ) • (l-Cx-Di) 



(9) 



: y — /£ 0 co 2 a em G r ef — icocc ee Vx G re f J • \1 — £o(o 2 a mm G re f — coa em Vx G re j 
: ^£ 0 co 2 a mm G re f + coa em Vx G re f) • (l — £ 0 co 2 a mm G re f — coa em Vx G re ^j 
- [in Q (i) 2 d em G K i + Kwa mm Vx G re f^ • (l — /i 0 co 2 a ee G re f — coa em Vx G re f^ 
= ^ 0 co 2 a ee G r ef + ffla em VxG re fj • (l — ii 0 a) 2 a ee G re f — coa em VxG re f^ 



(10) 

The Green's function can be written as an integral in the momentum space as 42 
1 

. ( n ) 



/oo 



Here, k 0 = k x x + k y y + k 0z z, k' 0 = k x x + k y y - k§ z z, e(k 0z ) =e(- k Qz ) = 
(xk y -ykx)/k//, h(k 0z ) = - (xk x + yk y )k Qz / (k Q k / /) + k//z/k 0 , h(-k 0z ) = 
(xk x +yky)k 0z /(k 0 k//) + k//z/k 0 and k/j = ^jk 2 x + k 2 y . The reflection coefficients 

Rj E and R™ are denned in the main text. 

In the calculations, G re f and its curl are first numerically evaluated and then one 
can apply equations (8-10) to evaluate the reflected fields E ref and H ref . Once the 
fields acting on the particle are obtained, one can then apply the analytic results in 
equation (3) to calculate the radiation pressure and spin density forces which are 
the source of the lateral force. 



Numerical simulation. All the full- wave EM simulations are performed with the 
package COMSOL Multiphysics. The relative permittivity of gold is described by a 



Drude model: £a u = 1 — ^p/^ 2 + K»co t ), where co p = 1.37 x 10 16 rads~ 1 
and co t = 4.084 x 10 13 rads~ 1 (ref. 60). For the simulation of the gold helix 
systems, we set / x w x t = 4 um x 4 urn x 200 nm and for the case of the dielectric 
substrate we set / x w x t = 4 um x 2.4 um x 200 nm. We choose a relatively 
larger substrate in the former case to reduce the effect caused by the reflection (at 
the edge) of the surface plasmons. For the simulation of the chiral sphere system, 
we set / x w x t = 2 um x 2 um x 200 nm. The finite substrate does not affect the 
physics here; the phenomenon also exists in the case of an infinite substrate. 
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